
#pragma once

#include <cmath>

//Function p
double p(double x);

//Function q
double q(double x);

//Function f
double f(double x);

//Assume that exact solution is available
double exactSolution(double x);

//Boudary Domain
const double a = 0; //vd1
const double b = 1; //vd1

//Value in 3rd Boudary condition
const double sigma0 = 0;    //vd1
const double sigma1 = 1;    //vd1

//Degree of Refinement
//The number of FE is 2^(MaxRefinemt)
//This value is needed when creating grid of FE
const int MaxRefinement = 3;

//Degree of aproximation of Gauss Integral
const int GaussDegree = 10;

const double PI = 3.1415926535897932384626433;
